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ABSTRACT 

It is shown that the diluted thermodynamical limit of a self-gravitating 
system proposed by de Vega and Sanchez suffers from the same problems as 
the usual thermodynamical limit and leads to divergent thermodynamical 
functions. This question is also discussed from the point of view of mean 
field theory. 
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Recently, de Vega and Sanchez have pointed out that a kind of thermo- 
dynamical limit of a self-gravitating system can be defined if one considers 
what they call the diluted limit: send the number of particles, N, and the vol- 
ume, V, to infinity, keeping constant the ratio N/V 1 ^ instead of the density, 
N/V [1]. This is a rather surprising -and very interesting- suggestion, since 
it is well known that the usual thermodynamical limit of self-gravitating sys- 
tems leads to singular thermodynamical functions (due to the gravitational 
instability) [2, 3]. However, we will show in this paper that the diluted 
regime leads to divergent thermodynamical functions in a way similar to 
the usual thermodynamical limit. Before demonstrating this statement, we 
will remember the differences, in what concerns the thermodynamical limit, 
between a thermodynamical stable system [4] and a self-gravitating system, 
and present the arguments that support the idea that the diluted limit gives 
well defined thermodynamical functions for self-gravitating systems. After- 
wards, we shall proof that this is not possible, and we will point out the 
loophole in the arguments that lead to the wrong conclusion. We will end 
the paper with a brief discussion on mean field theory. 

Consider a system of ./V classical particles, endowed with a hard core, 
the interactions of which are described by a potential energy which is a sum 
of two body contributions: 

^ N {r 1 ,...,r N ) = ^2</>(ri-rj). (1) 

i<j 

Using the notation fyj = <j)(ri — rj), the micro canonical and canonical par- 
tition functions, Zmc and Zc, respectively, can be written as 



N 

3AT/2-1 



JV i=l i<j 



<3 



(2) 
(3) 



where E is the energy, V the volume, j3 the inverse temperature, [x] + = x 
if x > and [x\ + = if x < 0, and we have ignored some irrelevant factors 
involving powers of the particle mass. Introducing the family of functions 

N 

G { y ] («) = / II ^ 6 ( u -J^H *V) . ( 4 ) 

JyN 1=1 i<j 
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the partition functions can be obviously written as 

^mc = mr ^ N/2) J duG^^iE-N-u}^ 2 - 1 , (5) 
J duG { y\u)exp[-f3N a u]. (6) 



Assuming that, as TV — ► oo, keeping the density p = N/V fixed, the 
following asymptotic behavior holds 

G { y\u) ^ e N9 ^P\ (7) 

we can get the partition functions with the aid of the saddle point method. 
For the microcanonical ensemble (MC) we have the equation 

dffO, P) = 3 1 (8) 
du 2 E/N -u K ' 

The solution is, obviously, a function of e = E/N and p = N/V, u = 
u m (e,p), and the entropy scales with N. The corresponding saddle point 
equation for the canonical ensemble (CE) is 

(0, 

the solution of which is a function of (5 and p, u = u c (f3,p), and the cor- 
responding thermodynamical potential (usually called the Helmholtz free 
energy) is extensive. It is not difficult to realize that the saddle point equa- 
tions imply the equivalence between the MC and the CE. 

Let us formally see that the scaling (7) holds if the two body potential 
is short ranged. We will use the notation 

</( ri ,..., rjv )) = / Hd 3 ri f(n,...,r N ). (10) 

^ i 

Using the Fourier representation of the Dirac delta we can write 

GP(u) = V N |^ e -GWH (11) 
where, with our definition (10), 

uj i r N uj 

G { y\u) = ( exp(-i^ hi) ) = / II ^ ex PH^ E hj) ■ 

i<j i=l i<j 

(12) 
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We can now apply the cumulant expansion to the above expression 



( exp(-iw ^2 fa) > = ex P [Yl Yl ■ ■ ■ Yl ■ ■ ■ &n*,)c 

i<j \n=0 h<jl in<jn 

' (13) 

where (4>i 1 j 1 • • • (fti n j n )c is the connected correlation function. It is well known 
that only those sequences of coordinates {hji, ■ ■ ■ ,i n jn} for which the in- 
tegral of Eq. (14) below cannot be split into the product of two (or more) 
integrals contribute to the connected correlation function. We will call such 
sequences connected sequences. For large N and n <C N, the number of 
connected sequences asymptotically is c n N n+1 , where c n is a number inde- 
pendent of N. Let us define £i = — rj t , for / = 1, . . . ,n. The number 
of connected sequences for which the £/'s are linearly dependent scales with 
a lower power of N, and, therefore, they do not contribute as N — ► oo. 
Only connected sequences with the &'s linearly independent contribute to 
the leading term in N of the n-th order of the cumulant expansion. For 
such sequences, if the potential decays sufficiently fast at long distances, 
the integral of (f> il j 1 ■ ■ ■ 4>i n j n over the &'s gives a number, Xhji-i„j n , that is 
independent of the volume if V is large. The integral over the remaining 
variables gives a power of the volume, V N ~ n . Therefore, we have: 



V 

yN J II " ' ' ■>»■><■ \ 

1=1 

The connected correlation function is a sum of products of integrals of the 
above type and has the same behavior. Collecting all the results, we have 

£ ■•■ E <«W ■ ■ *W»>c « AT* 1 ^ > (15) 

h<jl in<jn 

where Xn is the average of Xiiji-i n j n over the connected sequences. This 
behavior is broken down for n ~ N, but the contribution of these terms to 
the cumulant expansion is negligible. Then, the n-th term is 

CnXn (N\ n 

Performing the change of variables to — ► Nlo in the integral of Eq. (11), we 
get for a = 1 

Gy\u) « Yl %L e ^u + n^ p)] (17) 
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where 

oo 

n(x) = Y(-i) n ^x n . (is) 

n=0 

The integral of Eq. (17) can be evaluated by the saddle point method, and 
gives 

G^i^^^explNgiu/p)], (19) 

where 

g(u/p) = 0(u/p)u/p + n[oj(u/p)} , (20) 
and u)(u/p) is the solution of 

n'(w) = -«/p. (21) 

Therefore, if the potential is short ranged, G^(u) verifies the scaling of 
Eq. (7) and the thermodynamical functions are well behaved in the thermo- 
dynamical limit. 

What about a self-gravitating system? Since the potential decays as the 
inverse distance, we have that, in the same conditions as before 



N 

d Z T k Kh---4>i njn ^ Xi ^ njn , (22) 



where R is the linear size of the system (V = i? 3 ). The same arguments of 
the above paragraphs lead to the following expression for the n-th term of 
the cumulant expansion: 



(■ 



-')"S^^" /3 (f)"- < 23 » 



The above expression shows that, as expected, the asymptotic behavior (7) 
does not hold. However, we see that for a = 5/3 the n-th order term of the 
cumulant expansion is 

Hence, changing the variable u> — ► Nuj in the integral of Eq. (11) and using 
the saddle point method, leads to 

4 5/3) («) *^exp[A^(Vp 1/3 )] (25) 
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for a self-gravitating system. 

Now, the saddle point equation that gives the MC partition function is 

dgju/p 1 / 3 ) = 3 l 

du 2E/N^-u' { ' 

The solution is a function of E/N 5 / 3 and p, u = p 1 / 3 u m [E/(p 1 / 3 N 5 / 3 )], 
and therefore the entropy is not extensive 1 . The MC temperature, Tmc = 
{d\nZ MC /dE)-\ is 

TMC = W ~ l^P^E/iN^p^)] . (27) 

We see that the MC temperature diverges as the number of particles in- 
creases. This is a manifestation of the so called gravothermal catastrophe 
[5]. 

Let us analyze the canonical ensemble. If the temperature is fixed, the 
contribution of Gy (u) is negligible compared with the contribution of the 
Boltzmann weight, exp[— (3N 5 / 3 u], and the saddle point solution is given 
by u = u c = uo, where uq is the minimum possible value of u, that is, 
the minimum potential energy. Hence, the system is completely collapsed. 
The only way to avoid the complete collapse of the system as the number 
of particles increases is to increase the temperature by a factor iV 2 / 3 , in 
agreement with the MC analysis. 

The thermodynamical potentials for self-gravitating systems are non ex- 
tensive and the system will be inhomogeneous. The thermodynamical limit 
is singular and gives ill behaved thermodynamical functions. 

Up to here, everything is known. There is, however, another interest- 
ing possibility: as de Vega and Sanchez pointed out [1], it seems that the 
function Gy\u) will have the asymptotic behavior as exp[Ng(u, N/R)] if 
we keep a = N/R constant (instead of p = N/V) as N — > oo. They called 
this the diluted limit. Indeed, plugging V = R 3 in Eq. (23), the ra-th term 
of the cumulant expansion of Gy ^ (u) reads 



n\ N n ~ l \R / 

Then, Gy\u) m exp[Ng(u, a)], and the thermodynamical potentials will 
scale with N. In particular, the canonical partition function will scale as 

z c « Yl e -*m*) . (29) 



1 It is not an homogeneous function of N, E, and V. 
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The system will not be extensive, however, and will develop inhomogeneities, 
but the thermodynamical functions will be well behaved in the thermody- 
namical (diluted) limit [1]. 



The above proof of the existence of the diluted thermodynamical limit 
is formal, since it relies on a series expansion, the convergence of which has 
not been demonstrated. Indeed, we are going to show that the diluted limit 
cannot give well behaved thermodynamical functions. Let us consider the 
canonical partition function on a volume V = R 3 in the diluted regime, so 
that iV ~ R, and let us take a portion of such volume of linear size R$ < R, 

3 

o- 



such that N ~ Rq. Then, we obviously have 



N 



1 f V" 

Z C > w/- 3N/2 / N II^* exp[-/3 J>,-] > -J^exph/3 J>,> ] , 

Jv o 1=1 i<j ' i<j 

(30) 

where 

1 f N 

(^i)vb = yN N Yl d * rk hi ' ^ 31 ) 

and the last inequality follows from a well known property of the exponential 
function: (exp(y)) > exp((y)). Since {4>ij)v = — K /Ro, where k > is a 
geometrical number independent of Rq if Rq is large, we have 

V N 

Zc > h_ e PN(N^l)n/ Ra _ (32) 

Recalling that Rq ~ LqN 1 / 3 , where Lq is a fixed number with dimensions 
of length, we have 

N N 

Z C > -jfi- exp[N 5/3 PK/L \. (33) 

Therefore, the canonical partition function cannot behave as exp[— Nf(f3, a)] 
in the diluted thermodynamical limit. Rather, it behaves as in the usual 
thermodynamical limit. The reason is clear: the partition function is dom- 
inated by collapsed configurations even in the diluted regime. The gain 
in entropy provided by the dilution, which is of the order NlnV, cannot 
compete with the energy gain due to collapse, which of the order iV 5//3 . 

If the canonical partition function does not scale as exp[— Nf((3, a)], it is 
impossible that Gy\u) ~ exp[Ng(u, a)]. This is in conflict with the result 
of the cumulant expansion. The inequality (32) that originates this conflict 
is rigorous, so that the fallacy must be found in the cumulant expansion. 
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The solution of the paradox is that the cumulant expansion for G v {u) is 
dominated by terms of the order of N. It should converge for any finite N, 
but the radius of convergence shrinks to zero as N — > oo. Indeed, Gy\Nuj) 
is the analytical continuation of the Zc(P) to the imaginary axis. If the 
radius of convergence of the cumulant expansion were finite in the diluted 
thermodynamical limit, it would imply [cf. Eq. (17)] 

G { y ] (Nuj) «e™ M . (34) 

The asymptotic behavior of the above equation is not compatible with 
Eq. (33). The cumulant expansion for given by Eq. (18) must have 

a vanishing convergence radius in the diluted thermodynamical limit of a 
self-gravitating system. 

There is a clear explanation of the failure of the cumulant expansion. 
For imaginary u, G { y ] (Nuj) is a canonical partition function at temperature 
1/Im(a>). The cumulant expansion relies on connected correlation func- 
tions computed with a flat measure, V~ N FJ i d 3 ri, instead of the Boltzmann 
weight, Zq 1 Y\i d 3 ri exp[— Im(cj) ^ 4>ij\- The cumulant expansion, then, will 
be valid when the two measures are similar, i.e., in a gas phase. In a self- 
gravitating system, Eq. (33) indicates that the canonical partition function 
is completely dominated by the potential energy, and therefore that collapse 
takes place, at temperatures smaller or of the order of N 2 ^ 3 . Hence, for 
imaginary uj, the cumulant expansion will be only valid for Im(cj) < iV~ 2 / 3 . 
Thus, the convergence radius of the cumulant expansion shrinks to zero as 
iV~ 2 / 3 in the diluted thermodynamical limit. 

The collapse in the thermodynamical limit can be avoided by rescaling 
properly the radius of the hard core, a 2 . Indeed, it is known that in mean 
field theory the collapse phase covers the whole phase diagram in the limit 
in which the "filling" parameter, Na 3 /R 3 , tends to zero [6]. The "filling" 
parameter remains constant in the usual thermodynamical limit. To keep 
it constant in the diluted limit, the hard core must be scaled as a ~ N 2 / 3 . 
Then, previous argument does not apply, since such big particles do not fit 
in a volume of linear size R ~ N 1 / 3 . The minimum size of a region able 
to enclose the N particles must have a linear size scaling as N. Hence, the 
diluted thermodynamical limit can exist. This is, however, a rather trivial 
and ad hoc way of avoiding collapse: particles are forced to remain far away 
one from another, and it is physically difficult to justify the scaling of the 
hard core. Moreover, this way of preventing collapse in the thermodynami- 
cal limit cannot apply to other types of short distance regularizations, such 

2 I am grateful to P.-H. Chavanis and O. Fliegans for pointing this out to me. 
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as softened potentials [7]. Our proof of the non-existence of the diluted 
thermodynamical limit is also valid for softened potentials. 



It is generally believed that mean field theory is exact for self-gravitating 
systems. Since in mean field theory any thermodynamical function depends 
on the thermodynamical variables only through the dimensionless combina- 
tions A = -ER/(Gm?N 2 ) (MC) or 77 = (3Gm 2 N/R (CE), the thermody- 
namical limit must be taken keeping either A or rj finite. Then, it could be 
argued that the dependence of A or rj on N can be absorbed in any of the 
quantities entering A or ij. For instance, one could take R oc N and E oc N, 
as in the diluted limit. We will show in the following that this is not correct. 

Consider a softened gravitational potential, for instance —Gm 2 (j)(r)/R, 
with 4>{r) = R/Vr 2 + s 2 , where s is the softening scale. Dividing the finite 
region that encloses the system in W cells of linear size w (u> 3 = R 3 /W), and 
denoting by rij the occupation number of the ith. cell, the partition functions 
(2) and (3) can be approximated by [2]: 

ZMC = f(3iV72) £ rW [- A + 2>"^] + ' ( 35 ) 

^ ' > m=o n w =0 LU 1 i,j 



N N 



ZC = £, ^^exphyj^fe], (36) 

where 4>ij = 4>(ri — rj) and rj is the position of the center of the ith cell (we 
have again ignored constant factors). The above equalities hold rigorously 
in the limit W — > 00. For large N, the factorials can be approximated by 
their asymptotic form. Defining the density by pi = rii/(Nw s ), so that 
pi G [0, 1/w 3 ] becomes a continuum variable as N — > 00, we have, ignoring 
constant factors 

w 



Zmc = I \\ dp k 5(^2w 3 pi - 1) 

J /„ 1 7 



k=l 

3 

exp{iV[-u> 3 ^{Pi l ^Pi -Pi) + 2 ln [" A + Yl w& PiPM+]} > ( 37 ) 
W 



fc=i j 

exp{iV[-u> 3 ^(/^ln/?; - p^ + r/^u> 6 piPjfaj]} . (38) 



9 



Mean field theory is obtained by taking N — ► oo before W — > oo. In 
such case, the integrals in the above equations are saturated by the maxi- 
mum of the integrands, provided that the number of cells, W, and A or n, 
respectively, are kept constant. Obviously, the maximum of the integrand 
is given by the maximum of the entropy per particle, 

S = - J d 3 r[p(r)lnp(r) - p(r)] + | ln[-A + J d 3 rd 3 r' p{r)p{r')^{r - r')} , 

' " ' (39) 

in the MC, and with the minimum of the free energy per particle, 

T = - J d 3 r[p(r) lnp(r) - p(r)\ + 77 J d 3 rd 3 r' p(r)p(r')4>(r - r') , (40) 

in the CE, with the constraint f d 3 rp(r) = 1, if the grid is fine enough. The 
thermodynamical functions depend therefore on two dimensionless parame- 
ters: s/R and A or n. However, to derive mean field theory rigorously, one 
has to take the limit N — > 00 in (37) and (38), keeping the number of cells 
constant. Otherwise, if W grows with N, the integrals of Eqs. (37) and (38) 
cannot be evaluated by the saddle point. This implies that R must be kept 
constant. Therefore, the factors entering A and r] cannot be absorbed in 
R. They can be absorbed in G or in m? . Hence, mean field theory does not 
support the diluted limit either. 

In conclusion, the thermodynamical limit of a self-gravitating system 
does not exist, either in the usual form or in the diluted regime of Ref. [1]. 
The meaning of non-existence of the thermodynamical limit is that the ther- 
modynamical potentials do not scale properly with and thus thermody- 
namical functions, such as temperature, diverge. Nevertheless, it is possible 
to take the usual thermodynamical limit and, consequently, to use safely the 
usual thermodynamical tools by first regularizing the long distance behav- 
ior of the gravitational potential, introducing a very large screening length. 
The system is then thermodynamically stable and the thermodynamical 
limit does exist. Afterwards, one can study the limit in which the screening 
length tends to infinity [8]. 

I thank P.-H. Chavanis and O. Fliegans for some interesting discussions. 
This work has been carried out with the financial support of a Ramon y 
Cajal contract of Ministerio de Ciencia y Tecnologi'a of Spain. 
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